Depolarization for quantum channels with higher 
symmetries 



A B Klimov 1 and L L Sanchez-Soto 2 

1 Departamento de Fisica, Universidad de Guadalajara, 44420 Guadalajara, Jalisco, Mexico 

2 Departamento de Optica, Facultad de Fisica, Universidad Complutense, 28040 Madrid, 
Spain 



Abstract. The depolarization channel is usually modelled as a quantum operation that 
destroys all input information, replacing it by a completely chaotic state. For qubits this has a 
quite intuitive interpretation as a shrinking of the Bloch sphere. We propose a way to deal with 
depolarizing dynamics (in the Markov approximation) for systems with arbitrary symmetries. 
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Quantum systems evolve via unitary operations determined by the Schrodinger equation. 
This is also true for composite systems, e.g., a small quantum system S that is surrounded 
by some environment £, with which it interacts. However, while the evolution of the total 
system is described by a unitary operation Uss{t), the dynamics of the system S alone 
(obtained by tracing out the uncontrollable degrees of freedom of £) is, in general, no 
longer unitary (Breuer & Petruccione 2007, Weiss 2008, Alicki & Lendi 2007). The system- 
environment interaction leads to entanglement between them, and this is reflected into the 
fact that Use it) ^ Us (i) <8> Ue (t). From the perspective of quantum information processing, 
such an interaction is undesirable and causes errors and noise in the system (Nielsen & 
Chuang 2000). 

The dynamics of S can be described by a finite-time trace-preserving completely positive 
map (CPM) (Stinespring 1955, Choi 1972) that transforms input states g in = g(0) into output 
states £? out = g(t), i.e., 



which is also known as a quantum channel. One may think of the environment as extracting 
information from the system, as it will typically map pure states into mixed states. 

This noise process can be also described by a quantum operation involving only operators 
of the system of interest. This is called a Rraus decomposition (Kraus 1983) and has the form 
(omitting all the unnecessary subscripts) 



Qin I-* Qont 



£t(f?in) , 
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where the Kraus operators satisfy the condition 




(3) 



r 



which ensures that unit trace is preserved for all times. 
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A channel is Markovian when the coupling coupling of the system S with the 
environment £ can be treated under the Markov and Born approximations. The evolution of 
the state at a given instant is then fully determined by the state at that instant, so the process ha 
no "memory" of its past. This is a commonly used approximation in quantum optics and leads 
to the well-known Lindblad form of a master equation (Lindblad 1976, Gorini et al. 1976). 
For these Markovian channels, one can always write 

£ t (g) =e £t e(0), (4) 

where the Lindblad superoperator £ is 

C(g) = -i[H, g] + \ a - + [Uq, U}) • (5) 

r 

Here H is the Hamiltonian of the undamped system S, L r are system operators defined to 
model the effective dissipative interaction with the environment, and a r > are constants 
that account for decoherence rates. In this form, C appears as the generator of a CMP. In fact, 
expanding equation (5) to first order in the short-time interval r one can immediately find the 
corresponding Kraus operators (Shabani & Lidar 2005) 




(6) 

K r = \fr L r . 

For obvious reasons, the characterization and classification of these maps have attracted 
a lot of interest in recent years (Keyl 2002). The majority of the results obtained so far relate 
to two specific classes: the qubit channels and the bosonic Gaussian channels (Caruso & 
Giovannetti 2007). In this paper, we shall be mainly concerned with the former. 

In classical computation, the only error that can occur is the bit flip «-> 1. In quantum 
computation, however, the existence of superposition states brings also the possibility of other 
basic errors for a single qubit. They are the phase flip and the bit-phase flip. The first changes 
the phase of the state, and the latter combines phase and bit flips. The set of Kraus operators 
for each one of these channels is given by 

k = y/l-p/21, Ki = y/j/2aj, (7) 

where j = x gives the bit flip, j = z the phase flip and j = y the phase-bit flip. They can 
also interpreted as corresponding to a probability 1 — p/2 of remaining in the same state, and 
a probability p/2 of having an error. It is not difficult to represent these channels in terms of a 
master equation of the form (5): the associated Lindblad superoperator L{ turns out to be the 
corresponding Pauli matrix &j. 

The environment-induced noise is modelled using various simplified approaches. Here, 
we concentrate in the depolarization channel, which employs unbiased noise generating bit 
flip errors and phase flip errors and is represented by the CPM (Bennett et al. 1997) 

P 

£{0) = + o z ga z ) . (8) 

Note, however, that, since the three interaction channels corresponding to bit error, flip error, 
and phase error do not commute, one could rightly argue that the incoherent addition of these 
channels in the CMP (8) is, at least, questionable. 
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The Kraus decomposition of this channel reads as 

k = v/I - pi, 

(9) 

v p ~ v p ~ v P - 

while the associated Lindblad equation is 

§ = -r(^g-^ij , (io) 

where we have omitted the free evolution of the system, since it is irrelevant for our purposes 
here, and T is a constant. If we use the standard Bloch parametrization for the density matrix 

g = I (i + «.*), (ii) 

we immediately get 

s = -Ts, (12) 

which clearly shows that the action of the channel is to contract the sphere with a lifetime 
r _1 . This is precisely the idea behind depolarization: for long times the system will end in a 
fully depolarized or chaotic state, whose density matrix is diagonal g unpo i = \ 1 • 

For n qubits, a possible extension is to assume that the error operators can be represented 
by the n-qubit Pauli group (Nielsen & Chuang 2000) 

V n = {i,a x ,a y ,a z }® n , (13) 

where ®n denotes the n-fold tensor product. This means that each qubit is acted by the 
identical independent depolarizing channels. However, alternative models, such as collective 
or correlated depolarization, have been proposed (Reina et al. 2002, Banaszek et al. 2004, Ball 
et al. 2004, Ball & Banaszek 2005, Ban & Shibata 2006), which do not fit in this simple 
extension. The point we wan to stress is that the final result depends obviously of the 
symmetry of the channel and this, in turn, depends on physical considerations. 

The purpose of this work is to address this problem and identify the proper depolarizing 
channel for more general situations. To this end, consider a system whose Hamiltonian is 
a function of a symmetry Lie algebra 21. In order to maintain the discussion as simple as 
possible, we assume that this algebra is semisimple and denote by A the set of nonzero 
roots (Erdmann & Wildon 2006). We use the standard Cartan-Weyl basis {hi, e a } in terms of 
which we have the commutation relations 

[hi,hj] = 0, [hi,e a ]=a(hi)e a , [e a , ip] = N af3 e a+f} , (14) 

where the last one is valid only when a+[3 e A. The operators {hi} (i runs from 1 to £, where 
t is the rank of the group) constitute the Cartan subalgebra and may be taken diagonal in any 
irreducible representation. On the other hand, {e a , e_ Q } are raising and lowering operators 
and we can always choose = e_ a . 

The Hilbert space decomposes also into finite-dimensional subspaces 

H = 0^a (15) 

A 

so that in each H\ the operators of 21 act irreducibly. Let |h; A) be an orthonormal basis in 
H\, where h = (hi, . . . , hj, . . . , hi). Then, we have that 

hj\h; A) = hj\h; A) . (16) 
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Notice also that the density operator of any quantum unpolarized state can be written in terms 
of these invariant subspaces 

£unpoi = (D^a 1a , (17) 

A 

where 1 \ denotes the unity in the corresponding subspace and all the coefficients r\ are real 
and nonnegative and they fulfill the unit-trace condition. 

Since all the error operators are in the algebra 21, the only possible Kraus operators are 
also elements of 21, and they induce a local Lindblad equation fully analogous to (5). Very 
different processes can be described in terms of this master equation. The first one is what 
can be called a pure dephasing channel, represented by 

L r = h r . (18) 

Such a map obviously preserves the diagonal operators: C{h r ) = 0, and the asymptotic limit 
of the channel is given by 

lim £ n (g) = yWA)|h;A><h;A|, (19) 

n— >oo z — ' 
h,A 

where £»hh(A) are just precisely the occupation probabilities in each invariant subspace. For 
the case a symmetry algebra su(2), generated by {S x ,S y , S z }, the archetypal example of this 
situation is 

C(q) = S z gS z - q + e S 2 z ) = S z gS z - q, (20) 

which is a standard way of accounting for process that lead to a loss of coherence without 
changing the level populations (Briegel & Englert 1993). 

Next we consider a generalized amplitude-damping channel, for which 

L r = &-a T ■ (21) 

In the asymptotic limit, this leads to the ground state in each invariant subspace 

lim £ n (g) = V |h min ; A)(h min ; A| , (22) 

n^oo z — ' 
A 

where |h m ; n ; A) is the lowest-weight state that is annihilated by all the lowering operators: 
e- a |h m i n ; A) = 0. It is clear that now no Lindblad (neither Kraus) map preserves the diagonal 
operators. In particular, for su(2) symmetry, the textbook example is 

£{q) = S-qS+-±(S+S-q+qS + S-). (23) 

which is a standard model for damping (Agarwal 1974). 

While these generalizations are more or less obvious, the corresponding one for a 
depolarization channel is far from trivial. We claim that such a situation must be described 
by the action of a Lindblad operator proportional to e_ Qr followed by other proportional to 
e +ar . One can check that only in this way we asymptotically get an unpolarized state 

lim £ n {g) = V Tr[g(X)] t x = g unpoh (24) 

A 

where the trace operation is taken in each invariant subspace. The associated Lindblad 
operator is 

£($) = ~ XX ^-otrQZar + 2e Qr ge_ ar - {e Qr , e_ Qr } g - g {e ar , e_ Q? .}) , (25) 



Depolarization for quantum channels with higher symmetries 



5 



and one can check that C(g unpo i) = 0. An example of this situation for qubit systems parallels 
completely (25), but with the corresponding spin-like operators is 

C(q) = §-Q§+ + §+q§--^ {{§+,§-}§ +g{§+,e-}) . (26) 

It is worth noting that such a Linblad operator appears as a limit case of decaying into a bath 
at infinite temperature. 

An effective depolarization channel is actually a common situation in driven dissipative 
systems. Consider a typical evolution equation 

§ = -ig[§ x ,Q]+l[§-Q§+-± ({§+,§-} q+q{§+,S-}) , (27) 

which describes the decay (with rate 7) of an externally-driven (with coupling constant k) 
collective spin into a zero-temperature bath (Agarwal 1974). In the strong-pumping limit 
(g 3> 7), the above equation can be easily diagonalized (Klimov & Sanchez-Soto 2000): it is 
enough to apply the rotation U = exp(iirS y /2) and to make the rotating wave approximation. 
The final result is 

hd= ~ i9&, Qd] + \ (2S z qS z - S 2 J - gSfj 

+ 1(2S-qS + + 2S + qS--{S + ,S-}q-q{S+,S-}) , (28) 

where Qd = U g W is the density matrix in the rotated frame. We can clearly observe the 
emergence a pure dephasing (first line) and a depolarizing channel (second line). 

In summary, what we expect to have accomplished in this paper is to provide a 
construction of the depolarizing channel for systems with arbitrary symmetries. This may 
be more than an academic curiosity for more involved systems currently under investigation 
as candidates for quantum information processing. 
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